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The unimpeded relativistic propagation of cosmological nentrinos prior to recombination of the 
baryon-photon plasma alters gravitational potentials and therefore the details of the time-dependent 
gravitational driving of aconstic oscillations. We report here a first detection of the resulting shifts in 
the temporal phase of the oscillations, which we infer from their signature in the Cosmic Microwave 
Background (CMB) temperatnre power spectrnm. 
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Introduction: The hot, dense conditions of the early 
universe included a thermal background of photons. Ini¬ 
tially scattering with a high rate off of the free electrons 
in the plasma, they eventually streamed freely as the 
plasma cooled to ksT ~ 0.3 eV, the density of free elec¬ 
trons plummeted, and the mean free path became larger 
than the extent of the observable universe. Today we ob¬ 
serve these photons, about 90% of which last scattered at 
this epoch, as the cosmic microwave background (CMB), 
stretched by expansion into millimeter wavelengths. 

These same hot and dense conditions led to a cosmic 
neutrino background (CNB) that contributes nearly as 
much as photons to the total energy density in the early 
universe. The neutrinos began to stream freely at kgT ~ 
MeV, and continue to stream freely through the cosmos 
to this day. Unlike with photons, direct detection of the 
CNB is exceedingly difficult [1]. 

We have had indirect, though highly significant, evi¬ 
dence for the CNB for decades, starting with determina¬ 
tions of the primordial abundance of Helium, an abun¬ 
dance affected by the contribution of neutrinos to the ex¬ 
pansion rate (see Steigman [2]). More recently we have 
evidence from the amplitude of the damping tail region 
of the CMB angular power spectrum 0 a 0. 

Due to the indirect nature of these inferences of the 
CNB, they are highly model dependent. Therefore we 
are motivated to find signatures of the CNB that are es¬ 
pecially robust to changes in model assumptions. We 
report here the first detection of a particularly robust, 
though subtle, effect of the CNB: its influence on the 
temporal phase of the acoustic oscillations of the primor¬ 
dial plasma. 

Bashinsky and Seljak [B] analytically found that the 
propagation of neutrino perturbations at speeds faster 
than the speed of sound in the plasma alters the time- 
dependent gravitational driving of the acoustic oscilla¬ 
tions. Altering the driving of a harmonic oscillator has 
two effects: it changes the amplitude and the temporal 
phase of the oscillations. To date, no one has isolated 
these two distinct effects. Although the influence of neu¬ 
trino perturbations has been detected mm , it remains an 
open question whether the phase shift plays a significant 


role in these determinations. 

Isolating the phase shift effect clarifies the amount of 
model dependence in our inferences of the CNB. The am¬ 
plitude effect has a very smooth dependence on angu¬ 
lar scale, and thus can be mimicked by a smooth scale- 
dependent alteration of the statistical properties of the 
initial conditions. In contrast, mimicking the phase shift 
with an alteration to the initial conditions would require 
an oscillatory scale dependence that just happens to have 
a period matched to that of the acoustic peaks. It is also 
very difficult to mimic the phase shift by alterations to 
the matter content [6], short of substituting some other 
dark and relativistic species in place of the expected cos¬ 
mological neutrinos. 

In what follows we use a phenomenological parameter¬ 
ization of the influence of the phase shift on observables 
to directly determine the size of the phase shift due to the 
CNB. Although it is a subtle effect, with the Planck data 
it is finally possible to unambigously isolate this robust 
signature of neutrino free streaming. 

The power spectrum is Ci = (|a^rnP) where T{9,(p) = 
is a decomposition of the temperature 
map into spherical harmonics and (...) indicates an en¬ 
semble average. A undergoes I oscillations every 
360 degrees, so, e.g., Fisom are patterns with hot and 
cold spots separated by about 1°. We show model and 
measured power spectra in Fig. 1. 

The series of peaks is what we expect from a collec¬ 
tion of standing waves, all beginning their amplitude 
oscillations with zero initial momentum[14j . For exam¬ 
ple, the Fourier mode with spatial dependence 5pj oc 
cos(fc • x) will have a time dependent amplitude approx¬ 
imately given by Acos[fcrs(t)], where rs{t) is the dis¬ 
tance a sound wave can travel by time t (and we have 
ignored, for now, a time-dependent driving effect that 
alters this simple-harmonic motion solution). Modes 
with (comoving) wavelengths of ~ 400 Mpc project into 
i kd = 4 oo^Mpc ^ “ where d is the (comov¬ 
ing angular-diameter) distance to the last-scattering sur¬ 
face. These modes oscillate with a period such that they 
achieve their first extremum at the epoch of last scat¬ 
tering, t = U, and thus produce the first peak in the 
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FIG. 1: Undamped power spectra ICi (defined in Eq. with diflerent values of N^. In all panels the baryon density uib, 
the ratio of matter to radiation density pmlip^y + Pu) and the angular size of the sound horizon 9a are held fixed as these 
are well-determined by CMB data fairly independently of the assumed value of N^. In the top panel the dominant source 
of variation is the change in the damping scale 9d caused by the changes in N^. In the middle panel we fix 9d by varying 
the primordial fraction of baryonic mass in Helium appropriately, leaving the dominant source of power spectrum variation as 
the change in oscillation amplitude A'. Finally in the bottom panel, with the spectra normalized to remove the effect of A' 
variation, one can see the subtle impact of the shifts in temporal phase (f). The data points are the 2013 Planck data. 


power spectrum. Modes with a wavelength of ~ 220 
Mpc project into £ ~ 400. They have a shorter oscilla¬ 
tion period so that by t = t* they have gone through their 
first extremum and reached a null at the time of decou¬ 
pling. Modes that hit their extremum at decoupling 
contribute to the p**' peak. 

A key length scale for understanding the response of 
Cf, to the CNB is the sound horizon at decoupling, rs{t*). 
The sound horizon is smaller than it would be without 

the presence of the CNB, because the fractional expan- 
- 1/2 

sion rate H oc . If we scaled up the fractional expan¬ 
sion rate H at all times by a factor a, the decrease in 
time it takes for the temperature to drop to T ~ 0.3 
eV would lead to oc 1/a. The angle subtends, 

= 'rs{t^)/d where d is the angular-diameter distance 
to the last-scattering surface, strongly influences the lo¬ 
cations of the acoustic peaks such that 5£p = £pS9s/0s- If 
we knew d, we could use this effect alone to measure the 
energy density in the CNB. However, d depends on the 
(otherwise unknown) value of the cosmological constant. 

In Fig. we show a series of plots where we vary N^, 
while holding certain other quantities fixed, in order to 
demonstrate the observable consequences of various ef¬ 
fects of neutrinos. Because 0s, baryon density ujb, and 
the ratio of matter to radiation density Pm/{p~f + Pi^) 
are well determined by the data, in all rows we show 
variations with these parameters fixed. In the top row 


one can see the impact of N^, on the typical distance a 
photon diffuses prior to last scattering, rd- This diffu¬ 
sion suppresses anisotropy for modes with wavelengths 
A < r£), with an approximate effect of —>■ DiCi where 
Dg ~ exp [—where Ogj = ru/d. Because the 
diffusion is a stochastic process it scales with expansion 
rate as \fot rather than 1/a as Ts does. These different 
scalings mean that while we adjust d to keep 0^ fixed, we 
get 9d o: \foi. Thus for N^, = 5, the expansion rate is 
greater, leading to larger 9y) and more damping. 

To visualize more subtle effects of the CNB, we can 
vary the priordial fraction of baryonic mass in Helium, 
T/, to keep 0 d fixed as well [6]. Doing so in the mid¬ 
dle panel, we can see an impact of the perturbations in 
the CNB. As an initially over-dense region compresses 
under the influence of gravity, the compression does not 
occur rapidly enough to prevent the gravitational poten¬ 
tial from decaying due to the expansion-driven drop in 
density. By the time pressure gradients halt the compres¬ 
sion, the gravitational potential has nearly completely de¬ 
cayed. The result of this temporary time-dependent grav¬ 
itational driving of the accoustic oscillations is a change 
to the subsequent amplitude and phase so that the ampli¬ 
tude of a standing wave is A! cos[krs{t) + />]. The values 
of A! jA and depend on the details of the potential de¬ 
cay, and in particular on the fraction of the radiation 
that can freely stream out of over densities at the speed 
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of light. We can see in the middle row that increasing 
Nij decreases A'/A. 

Finally, in the bottom row we normalize the spectra to 
remove the effect of changing A'/A so that the change in 
values of (j) is more evident. As (j) changes, the value of 
k = kp for which kprs{tAj + (j) = pn changes so £p ~ kpd 
changes. The net result is Sip = —Scji/Os. In this Letter 
we show that these subtle shifts are detectable with the 
Planck data. 

Template fitting: To quantify the sensitivity of the data 
to the expected phase shift we must be able to artificially 
increase or decrease it, independent of other effects of 
the CNB. Since the phase shift effects are most observ¬ 
able deep in the damping tail, we work with an (approx¬ 
imately) undampted spectrum 

lCe = l{i+l)/{2TT)CiDf\ (1) 

where approximately undoes the damping expected 
for three neutrino species and standard Big Bang Nucle¬ 
osynthesis. We then define a transformation that sends 
1C I —> ICi+se^y which is controlled by a new parameter 
Following the analytic work of Bashinsky and Sel- 
jak m we take the amplitude of Sii, to be linear in the 
fraction of radiation energy density in neutrinos, R{N^), 
so that 

6i, = A{Nt*,N,)fu ( 2 ) 

with defined such that 

A{Nt^, N,) oc R{Nt^) - R{N,), (3) 

with proportionality chosen to match the shift from the 
standard 3.046 to one neutrino specie. 

The analytic result of [5] was obtained by ignoring the 
presence of any pressureless matter, an approximation 
that becomes increasingly valid as one considers modes 
that began oscillating deeper in the radiation-dominated 
era. With this approximation they find a Sf (and there¬ 
fore /f in our language) that is independent of k and 
therefore independent of i. To avoid this approximation 
we have numerically calculated Sii, from the Ci calcu¬ 
lated with the Boltzmann equation solver CLASS |H] for 
a suite of models with varying Ni,. To define the suite 
of models we compare two cosmologies: one from the 
ACDM posterior of the 2013 Planck data, and the other 
with nonstandard N^, drawn from 1 < N^, < 6. We iso¬ 
late the effects of neutrinos from the effects of the other 
components contributing to gravitational potentials, as 
well as effects of changing N^, to the background Hubble 
rate H{z), by following Hou et al. [3] in hxing the baryon 
density Wf,, redshift of matter-radiation equality Zeq, Silk 
damping scale On, and sound horizon scale 0s, as well as 
the spectral amplitude and tilt describing the adiabatic 
initial conditions from inflation between the two models. 
Since contributions to anisotropy arising after recombi¬ 
nation project differently, we zero out contributions from 
the integrated Sachs Wolfe (ISW) effect. 



FIG. 2: Posterior reconstructions of the phase shift template 
fi = a ln(Z) -I- d under the ACDM -I- Ai, -|- model, nor¬ 
malized to give the amplitude of the relative shift between the 
fiducial 3.046 neutrino species and a single specie, along with 
numerically obtained phase shifts obtained from 100 different 
simulated cosmologies with varying neutrino species drawn 
from Ni, = [1,6] relative to the fiducial cosmology, rescaled 
by the amplitude parameter given in equation Also in¬ 
cluded is the analytic approximation given by [B] , which they 
also found to be ~ 25% higher than the numerical result at 
t ~ 3000. 

We sample 100 different cosmology pairs, and find that 
in the region of parameter space explored by these mod¬ 
els, the phase shift is well captured by a linear response 
proportional to the fraction of radiation density in free- 
streaming neutrinos. The multipole dependence is well 
described by a logarithmic template, which is jointly sam¬ 
pled with cosmology against both the March 2013 Planck 
temperature data and the measured phase shifts in the 
100 cosmology pairs. Posterior samples of the template 
are shown in Fig.[^ along with the measured phase shifts 
Slv at the peaks in the ISW-less temperature power spec¬ 
trum JCf"^, as well as, for visualization, the locations 
of the corresponding peaks in the polarization spectrum 
ICf^, which is not used in the fit. 

To apply this template to models we confront with 
the data, we first decompose the temperature power 
spectrum C; into ISW and ISW-less components: C/ = 

(^noISW (Jisw (JCTOSS^ (jcross ^ (nO 

ISW) contribution. The artihcial space shift we intro¬ 
duce is given by 

ICi -)■ (4) 

5i,. = A{Nt^K)f{L), 

with the implied definitions for the various ICc compo¬ 
nents. The artificial £—space shift alters the power spec¬ 
trum away from that of the physical model which always 
has Ny = Nff, for which A = 0. To vary just the phase 
shift effect, we can set W = 3.046, the fiducial value, 
while varying N^'i. 
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FIG. 3: Top: 2D constraints on the jointly varying 

ACDM+A^^+A^^"^ parameter space. The constraints on 
(damping) and (phase shift) are essentially orthogonal. 
Bottom: Constraints from March 2013 Planck temperature 
power spectrum measurements on the number of neutrino 
species from (1) blue/solid: varying while holding Ni^ 
fixed at three and (2) green/dashed: varying along the physi¬ 
cal direction Ni, = . The constraints assume a Gaussian 

r prior of mean — 0.085 and width a = 0.015. 


Results from Planck: We use the publically available 
likelihood code clik [3] to determine constraints from 
the 2013 Planck temperature power spectrum measure¬ 
ments, with the polarization constraints approximated as 
a Gaussian prior on the optical depth to last scattering 
T for simplicity [TS]. We place uniform (flat) priors on 
the parameters N^, and which results in a flat prior 
for the physical case where N^, = We explore the 

model space using the MCMC routines provided by the 
Python library CosmoSlik m- 

In Fig. we show the constraint on from the 
ACDM model. For comparison, we include the 

constraints on the ACDM +N^ model space, with the 
phenomenological amplitude , N^) = 0. We find 

best-fit values of N^, = 3.3^q 2 [E] and 
To quantify the signihcance of the detection we con¬ 
ducted a likelihood ratio test. We found that the min¬ 
imum {= —21nL) decreased by 19.9 when switching 
from a ACDM model with = 0 to the ACDM -f 
model. A difference this large or larger will occur, as¬ 
suming = 0, with the same probability, 8 x 10 as 
a 4.5(7 Gaussian fluctuation. 

While letting AA and vary independently (top of 


Fig.i , the width of the constraint on Ni, is dominated 
by an rig- degeneracy: fixing Ug results in roughly a 
halving of the characteristic width of the posterior in the 
Ni, direction. For the direction, no such correlations 
exist, so while our constraints on depend somewhat 
on the characterization of initial conditions due to infla¬ 
tion (as well as, in extended models such one with the 
helium fraction Yp free), our constraint on is due to 
a feature in the data that is difficult to mimic with other 
cosmological degrees of freedom. In addition, we find the 
constraints on N^, and are nearly uncorrelated. This 
lack of correlation follows from the fact that the response 
of Cl to changing the phase shift is essentially orthogo¬ 
nal to the response due to other observable effects of the 
cosmic neutrino background in the CMB. 

Finally, we note that there is a slight dependence on 
priors for the 2D posterior shown in Fig. If we switch 
from uniform priors on the number of neutrino species AA 
and to their corresponding neutrino fractions R{N^) 
and R{Nl’f‘), the average value of the posterior shifts 
down by (AA^j, = 0.3, = 0.5), a shift of slightly 

more than 0.5 (t in both directions. This is predominantly 
due to a contraction in the high probability region at high 
Ni, or Regardless of prior, = 0 is heavily dis¬ 
favored. 

Conclusions: In this letter, we present the first detec¬ 
tion of the temporal phase shift generated by neutrino 
perturbations during the acoustic oscillation phase of cos¬ 
mological evolution, and find an amplitude of this effect 
consistent with the standard value associated with the 
three known neutrino species. As pointed out by Bashin- 
sky and Seljak [5] this subtle signature is a particularly 
robust one, difficult to mimic by either changes to the ini¬ 
tial conditions or matter content. Our detection of this 
effect is thus the most model-independent determination 
to date of the existence of the CNB. CMB Polarization 
measurements will bring even more robustness as they 
are capable of definitively ruling out the remote possi¬ 
bility of mimicking the phase shift effect by a designer 
alteration of the initial conditions. 

Were we to have found Ni, ^ we would be com¬ 
pelled to look for a physical explanation, such as v n 
scattering which could inhibit free streaming mm- 
With current data, we see consistency with the standard 
model. Future datasets, such as from Planck including 
polarization spectra m and CMB stage-4 power spectra 
m will provide even stronger constraints on the phase 
shift, either providing a signature of new physics or in¬ 
creasing conhdence in the standard cosmological model. 
We forecast that these will achieve = 0.41 and 

0.09 respectively. 
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